Highly-charged ions as a basis of optical atomic clockwork of exceptional accuracy 
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We propose a novel class of atomic clocks based on highly charged ions. We consider highly- 
forbidden laser- accessible transitions within the 4/ 12 ground-state configurations of highly charged 
ions. Our evaluation of systematic effects demonstrates that these transitions may be used for 
building exceptionally accurate atomic clocks which may compete in accuracy with recently proposed 
nuclear clock. 
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Atomic clocks are arguably the most precise scientific 
instruments ever built. Their exquisite precision has en- 
abled both foundational tests of modern physics, e.g., 
probing hypothetical drift of fundamental constants 
and practical applications, such as the global position- 
ing system. State of the art clocks carry out frequency 
measurements at the eighteenth decimal place [2j. As 
the projected fractional accuracy of such clocks is at the 
level of 10~ 18 Q, 0], it is natural to wonder how to ex- 
tend the accuracy frontier even further. We are only 
aware of one proposal, the nuclear clock that holds 
the promise of reaching the improved 10~ 19 accuracy 
level. The nuclear clock, however, relies on a yet unob- 
served optical transition in the radioactive 229 Th nucleus. 
Here we show that the nuclear clock performances can be 
replicated with atomic systems, fully overcoming these 
challenges. We identify several highly-forbidden laser- 
accessible transitions in heavy stable isotopes of highly- 
charged ions (HCI) that may serve as clock transitions. 
Similarly to the singly-charged ions of modern clocks Q , 
HCIs can be trapped and cooled d, 0] • The key advan- 
tage of HCIs comes from their higher ionic charge. As 
the ionic charge increases, the electronic cloud shrinks 
thereby greatly reducing couplings to detrimental exter- 
nal perturbations. Our analysis of various systematic 
effects for several HCIs demonstrates the feasibility of 
attaining the 10 -19 accuracy mark with existing technol- 
ogy 

Atomic clocks operate by locking the frequency of 
external oscillator (e.g., laser cavity) to a quantum 
(atomic/nuclear/ionic/molecular) transition. One tells 
time by simply counting the number of oscillations at 
the source and multiplying it by the known oscillation 
period. A suitable clock transition should have a good 
quality factor (Q-factor). Moreover, the clock transition 
frequency must remain unaffected by external perturba- 
tions. Finally, one has to be able to interrogate quan- 
tum oscillators for a long time, so the atoms should be 
trapped. The clock stability and accuracy generally im- 
prove with higher frequency of the clock transition, I'ciock, 
and the current accuracy record [2] is held by singly- 
charged ion clocks operating at optical frequencies. 



Before we start with the clock estimates, we would like 
to recapitulate a few basic facts about HCIs. In a multi- 
electron atom, optical electrons move in mean-field po- 
tential created by other electrons and the nucleus. How- 
ever, as the electrons are stripped away from the atom, 
the field experienced by the optical electrons becomes 
increasingly coulombic, and one could invoke intuitive 
hydrogen- ion-like estimates For example, the size of 
the electron cloud scales with the residual nuclear (ionic) 
charge Zi as 1/Zi. Since typical matrix elements are 
proportional to some power of atomic radius, most of 
the couplings to the detrimental external perturbations 
scale down with increasing Zi. Higher-order responses, 
e. g., polarizabilities, are suppressed even further due to 
increasing energy intervals that scale as Zf. Such sup- 
pression of couplings to external perturbations is the key 
to improved accuracy in the proposed HCI clock. 

Trapping and cooling clock ions beneficially increases 
interrogation time and reduces Doppler shifts. HCIs can 
be loaded in ion traps @, 0], however, due to the Zf 
scaling of intervals most of HCIs lack low-energy elec- 
tric dipole transitions that can be used for direct laser 
cooling. As in the state-of-the-art optical ion clocks 0], 
to circumvent this limitation, one may choose to em- 
ploy sympathetic cooling. In this scheme, long-range 
elastic Coulomb collisions with continually laser-cooled 
Be + ions drive HCI temperature down to mK temper- 
atures. Heavy HCIs may be cotrapped with relatively 
light ions of low ionic charge, because equations of mo- 
tion in ion traps depend only on the ratio of ion charge 
to its mass, Zi/M. For example, Ref. [6| experimen- 
tally demonstrated sympathetic cooling of Xe 44+ with 
Bc + ions. The basic idea Q is to initially cool HCIs 
resistively and then load precooled HCIs into the Be + 
ion trap. At sufficiently low temperatures the rates of 
undesirable charge-exchange processes between two ionic 
species become negligible. Heavier cooling species like 
Mg+ can be also used Q to improve mass-matching and 
thereby the cooling efficiency. Additional advantage of 
co-trapping two ionic species comes from the possibility 
of carrying out quantum logic spectroscopic clock read- 
out [2j and initialization. 
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FIG. 1: Systematic trend of basic clock properties for the 
[Pd]4/ 12 isoelectronic sequence. Clock frequencies (lower 
panel) and the ratios of radiative line width to the clock fre- 
quency (upper panel) are shown as a function of ionic charge. 
Clock frequencies lie in laser-accessible near-infrared and op- 
tical domains, while the narrow transition width assures high 
quality factor. 



There are many possible choices of ions. Including 
all degrees of ionization of the first 112 elements of 
the Mendeleev (periodic) table leads to 6,216 poten- 
tial ions. We are interested in those ions which have 
closed highly-forbidden transitions in the optical fre- 
quency band. Some of the optical transitions in HCIs 
were identified in Refs. [13-Gj]. We analyzed several 
possibilities, and we find HCIs with the [Pd]4/ 12 ground- 
state electronic configuration to be especially promising 
for precision timekeeping. 

The [Pd]4/ 12 configuration is the ground state con- 
figuration for all ions starting from Re 17+ which have 
nuclear charge Z > 75 and ionic charge Zi = Z — 58. 
We computed properties of such ions using relativistic 
configuration interaction method described in 13, ITij ; 
details of the calculations will be presented elsewhere. 
According to the Hund's rules, in these HCIs the 4/ 12 
and 4/ 12 ^4 states are the ground and the first ex- 
cited states respectively. Our computed clock frequencies 
^ciock and ratios of radiative width 7 to v c \ oc k are shown in 
Fig. Q] The transition frequencies range from the near- 
infrared to the optical region and are laser-accessible. 
The clock states are exceptionally narrow, as the upper 



TABLE I: Clock-related properties of representative HCIs of 
the [Pd]4/ 12 isoelectronic sequence. A c iock is the wavelength 
of the clock transition, r is the lifetime of the upper 4/ 12 ^4 
clock level, and Q is the quality factor. Systematic clock 
shifts are governed by differential static electric-dipole polariz- 
ability Aa E1 (0), black-body coefficient /3bbr and quadrupole 
moments of the clock states Q e . Hyperfine structure of the 
clock levels is determined by the nuclear spin I and hyperfine 
structure constants A. Numbers in square brackets represent 
powers of 10. 
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clock state may decay only via highly-suppressed electric- 
quadrupole (E2) transition. The resulting lifetime of a 
few hours leads to the prerequisite high Q-factor of the 
clock transition, so that 7/(27ri' c i ock ) remains below the 
10 -19 accuracy goal. 

Clock-related properties of several representative HCIs 
are compiled in Table |U As an example, below we fo- 
cus on 209 Bi 25+ ion. It has the highest transition fre- 
quency and Q-factor among stable isotopes. Being the 
heaviest among such isotopes has additional advantages 
as its large mass suppresses systematic effects related to 
Doppler shifts (see below). 

209 Bi has the nuclear spin of I = 9/2 and, as sketched 
in Fig. [2j the electronic states are split into a multitude 
of hyperfine components. The clock transitions must be 
insensitive to external perturbations, such as magnetic 
and electric fields. Because of that we choose specific hy- 
perfine states and magnetic substates: \F — 17/2, Mp = 
±5/2) o |F = 13/2, Mp = ±9/2) for clock transition. 
Similar to the virtual transition technique demonstrated 
in Hg + clocks the HCI clock operates on two tran- 
sitions that have opposite g-factors. Averaging over the 
two transition frequencies eliminates the linear Zeeman 
shift, making clock insensitive to B-fields. While such 
a technique could be applied to multiple transitions, as 
shown below, we further required that our specific choice 
minimizes couplings to electric field gradients. 

Clock accuracy is affected by multiple systematic ef- 
fects: magnetic fields, electric fields, Doppler (motion- 
induced) effects, blackbody radiation (BBR), and grav- 
ity. We show that all these effects are suppressed at the 
desired 10 -19 fractional accuracy. 

We start with examining the BBR shifts; these arise 
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FIG. 3: Differential quadrupole moment AQ as a function 
of magnetic quantum number of the excited clock level for 
209 Bi 25+ . Different curves are labeled by values of the total 
angular momentum of the excited state F. The ground state 
hyperfine component remains fixed \F g = 13/2, M g = 9/2}. 
The minimial value of AQ is encircled; it is attained for the 
\F e = 17/2, Me = ±5/2) "magic" state. 



FIG. 2: Proposed virtual clock transitions in 209 Bi 25+ highly- 
charged ion. Averaging over the two indicated transition fre- 
quencies removes the first-order Zeeman shift. Specific choice 
of hyperfine components and magnetic sub-states minimizes 
shifts due to couplings to gradients of the trapping electric 
field. 



due to perturbations by the photon bath at room tem- 
perature. The fractional shift reads [15| 
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where Aa(0) is the differential static polarizability of the 
clock transition and T is the BBR temperature. For 
HCIs, polarizabilities are suppressed as 1/Zf and our cal- 
culation yields Aa(0) ~ —2.3 x 10~ 4 ciq (arj is the Bohr 
radius). This tiny Aa(0) translates into the fractional 
BBR shift at room temperature of just 5.8 x 10~ 21 . Sim- 
ilarly, differential polarizability determines sensitivity to 
stray electric fields: A^/i/ c i oc k = -Aa(0)f 2 / (2/ii/ c i oc k). 
Typical E-fields of 10 V/m lead to negligible 10" 284 -level 
shifts. Cooling lasers shining on the coolant ion will lead 
to AC Stark shifts of the HCI clock levels. Again com- 
pared to the singly-charged ion clocks these shifts will be 
strongly suppressed due to much smaller Aa and also 
because the HCI and the coolant ion are repelled by 
stronger Coulomb forces reducing the overlap of the cool- 
ing laser beam with the HCI. 

The clock ion is trapped in a non-uniform field; the 
gradient of this field couples to the quadrupole moment 
Q of the clock states 16|. The quadrupole shift (QS) of 



the clock transition is given by 
A^qs _ 1 



^clock 



2/^clock OZ 



(1) 



where AQ ee (Qo) e — (Qo)g is the difference in expecta- 
tion values of the zeroth component of the quadrupolar 
tensor for the upper and lower clock states. 

The relevant clock shifts for singly- charged ion clocks 
are sizable and considerable efforts has been devoted to 
mitigating this effect [l|, [lfjj]. While for the HCIs one 
expects 1/Zf suppression of Q-moments, we find that 
the relevant clock shifts can be still appreciable. Below 
we minimize the quadrupole shift by exploiting the rich- 
ness of the hyperfine structure of the [Pd]4/ 12 HCI clock 
states. Indeed, the Q-moments of various hyperfine sub- 
states \aJI; FMp) attached to electronic state \aJ) may 
be conveniently expressed as a product of Q-moment of 
the electronic state Q e (aJ) and a kinematic factor 

(aJI]FM F \Qo\aJI;FM F ) = 

(3M| - F (F + 1)) K (J, I, F) Q e (aJ) 

where Q e are listed in Table Hand the Mp- independent 
factor K (J, /, F) reads 
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As the gradient is fixed by trap parameters, we mini- 
mize the difference (Qo) e — (Qo) g by considering all pos- 
sible pairs of hyperfine sub-states allowed by the E2 se- 
lection rules. The search for "magic" transitions depends 
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on the ratio of electronic Q-moments. These ratios can 
be determined experimentally by measuring frequencies 
of several hyperfine transitions in a trapped ion. In Fig. [3] 
we illustrate such a search based on our computed val- 
ues of Q-moments for 209 Bi 25 +. We find that the min- 
imal value of AQ = — 5 x 10~ 6 |e|aQ is attained for the 
magnetic components \F e = 17/2, M e = ±5/2) - \F g — 
13/2, M g = ±9/2). These are the clock transitions indi- 
cated in Fig. [21 

The field gradient in Eq. (|XJ) is fixed by the trap; we 
adopt the value from the Be + /Al+ clock [17|: d£ z /dz w 
10 8 V/m 2 . Additional gradient on the clock HCI is ex- 
erted by the coolant ion. For typical ion separations of 
10 //m the resulting gradients are smaller than the in- 
dicated trapping field gradient. With such gradients, we 
find Auqs /v c \ock ~ 5 x 10 -19 , which can be substantially 
reduced further. Indeed, due to rotational symmetry ar- 
guments, the QS can be fully zeroed out by averaging 
clock measurement over three orthogonal directions of 
quantizing B-field jltij . The power of this technique has 
been experimentally demonstrated for the Hg + clock l], 
where the QS was reduced by a factor of 200. The averag- 
ing out was not exact due to technical alignment issues. 
Combination of the "magic" choice of clock states with 
the averaging technique leads to our projected QS uncer- 
tainty of Av QS /v cio(ik = 2.5 x 10~ 21 . 

Clock frequencies are affected by magnetic fields. The 
first-order Zeeman shift can be eliminated by averaging 
the measurements over two virtual clock transitions indi- 
cated in Fig. [2l The dominant source of Zeeman-related 
uncertainties comes from AC B-fields caused by currents 
at the RF trap frequencies in conductors near the ion. 
Ideally, such fields would vanish at the trap axis, but in 
practice Bac is always present due to geometric imper- 
fections [H- We adopt Bag =5x10~ 8 T measured in the 
Al + /Bc + trap [H as the typical value. The AC fields con- 
tribute to the second-order Zeeman shift. Calculations 
of the relevant differential magnetic-dipole polarizability 
Aa M1 are dominated by intermediate states of clock hy- 
perfine manifolds. For our choice of magnetic substates 
for 209 Bi 25+ ; we find Aa Mi _ _2 i x io 10 Hz/T 2 which 

translates into a fractional clock shift of 4 x 10~ 20 ; it is 
below the sought accuracy goal. 

Working with HCI requires relatively high vacuum at- 
tainable in cryogenic traps cooled with liquid helium. In 
the context of ion clocks, such traps were demonstrated 
for Hg + ion [![. The rate coefficient [l8| for charge- 
exchange collisions of heavy HCIs with residual He atoms 
scales as Zc k pa 0.5 x 10 _9 ZjCm 3 /s. If the HCI were 
to survive for an hour, the number density of He atoms 
would have to be limited by 2 x 10 4 cm~ 3 . 

The zero-point-energy motion of trapped ion has a pro- 
found effect on the clock accuracy via the effect of special 
relativity, time dilation [l9[ . The fractional effect of time 
dilation may be evaluated as the ratio of the ion kinetic 



energy K to its rest mass energy, 

^TD^clock = ~K/Mc 2 . (2) 

To estimate the effects of time-dilation, we adopt trap 
parameters from the ion clock of Ref. [3] based on a pair 
of Al + /Be + ions. Indeed, once the trapping fields are 
specified, ion motion is entirely characterized by the ra- 
tio of ionic charge to its mass, Z^/A, where A is the 
atomic weight. For all the enumerated highly-charged 
clock ions, this ratio is about 0.15 which is comparable 
to the Zi/ A-r&tio for Be + . Since the trapping parameters 
are similar, the value of kinetic energy remains roughly 
the same as in the Al + /Be + clock, while the enumerated 
HCIs are about 10 times heavier than Al. This mass dif- 
ference leads to suppression of the time-dilation effects 
with heavy HCIs, see Eq.([2j). In the demonstrated 27 A1 + 
clocks the uncertainties due to time-dilation are at the 
level of a few parts in 10~ 18 with the goal of reaching the 
10 -18 accuracy milestone. Due to the mass scaling argu- 
ment we anticipate that 10~ 19 is the plausible accuracy 
goal for the proposed HCI clocks. Notice that this lim- 
itation is not fundamental as it relates to the technical 
ability to control stray electric fields in the trap. 

The clocks are affected by the effects of general rel- 
ativity as well 1 19]. The fractional frequency difference 
between two clocks at differing heights on Earth's sur- 
face is Avq I t'ciock = g Ah/c 2 , where g is the gravitational 
acceleration and Ah is the difference in clock height po- 
sitioning. If the two identical clocks differ in height by 
1 mm, the clock would acquire a lxlO -19 fractional fre- 
quency shift. Such uncertainty would limit accuracy of 
time transfer. 

To summarize, we have shown that the highly- charged 
ions may serve as a basis of optical atomic clockwork 
at the 10 -19 fractional accuracy. Such accuracy results 
from the smallness of electronic cloud in ions and there- 
fore suppressed couplings to external perturbations. The 
10 -19 fractional accuracy is matched only by the pro- 
posed nuclear clock [5|; our proposed clock avoids com- 
plications of radioactivity and uncertainties in transition 
frequencies associated with the nuclear clock. 
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